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ABSTRACT 

A small strain theory, three-dimensional elastic-plastic finite element analysis was used to 
simulate fracture in thin sheet 2024-T3 aluminum alloy in the T-L orientation. Both 
straight and tunneled cracks were modeled. The tunneled crack front shapes as a function 
of applied stress were obtained from the fracture surface of tested specimens. The stable 
crack growth behavior was measured at the specirhen surface as a function of applied 
stress. The fracture simulation modeled the crack tunneling and extension as a function 
of applied stress. The results indicated that the global constraint factor, a g , initially 
dropped during stable crack growth. After peak applied stress was achieved, a g began to 
increase slightly. The effect of crack front shape on a g was small, but the crack front 
shape did greatly influence the local constraint and through-thickness crack-tip opening 
angle (CTOA) behavior. The surface values of CTOA for the tunneled crack front model 
agreed well with experimental measurements, showing the same initial decrease from 
high values during the initial 3mm of crack growth at the specimen's surface. At the 
same time, the interior CTOA values increased from low angles. After the initial stable 
tearing region, the CTOA was constant through the thickness. The three-dimensional 
analysis appears to confirm the potential of CTOA as a two-dimensional fracture 
criterion. 


NOMENCLATURE 

b Uncracked ligament length, mm 

B Sheet thickness, mm 

v Poisson's ratio 

E Elastic modulus, MPa 
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a Half crack length, mm 

S Applied stress, MPa 

Aa s Crack growth at the specimen's surface, mm 

x,y,z, Cartesian coordinates, mm 

u,v,w Displacements in the x, y, and z directions, mm 

On Normal stress in the i direction, MPa * 

Oq Flow stress (average of yield and ultimate), MPa 

a g Global constraint factor 

A[ Projected area of element i on the crack plane, mm 2 

A m Total projected area of all yielded elements on the crack plane, mm 2 

K Stress-intensity factor, MPa Vm 

K n Normalized stress-intensity factor, MPaVm 

\|/ c Critical crack-tip opening angle, degree 

Vj v-displacement at location of \jr c measurement, mm 

Xi Distance from crack-tip to location of y c measurement, mm 

w Half width of specimen, mm 

INTRODUCTION 

The prediction of residual strength in an aircraft fuselage structure requires a 
fracture criterion that can account for large amounts of stable tearing in thin-sheet 
material. One criterion that has been proposed is the crack-tip opening angle (CTOA) [1 - 
6]. This criterion states that stable crack growth will occur when the local angle made by 
the upper crack surface, crack-tip, and lower crack surface reaches a critical value. 
Experimental studies have shown that at the surface, initial CTOA values are high, but 
quickly decrease and maintain a constant value after a small amount of stable crack 
growth [7-9]. Reuter and Lloyd [10] have shown experimentally that for a surface crack, 
CTOA is not constant around the crack perimeter. Two-dimensional numerical studies 
that neglect the initial transient and through-thickness variation in CTOA behavior have 
shown that this criterion can successfully be used to predict residual strength for a 
number of different specimen configurations [8, 11-13]. 

The early stage of stable tearing where the measured surface CTOA is not 
constant has prompted concerns about the local crack-tip stress fields and the constraint 
in thin sheet materials. Horn and McMeeking [14] analyzed stationary, straight, through- 
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thickness notches in thin sheets using a large deformation theory, three-dimensional 
elastic-plastic finite element analysis. Newman et al [8, 15-16] analyzed stationary, 
straight and tunneled cracks, using a small strain theory, three-dimensional elastic-plastic 
finite element analysis. Both found that considerable constraint develops even in thin 
materials where significantly higher stresses develop in the interior than on the surface. 

In addition to constraint variations and the large CTOA values during early stages 
of crack growth, significant changes in crack front shape also occur. Sutton, Dawicke, 
and Newman [17] have shown that for 2024-T3 aluminum alloy in the T-L orientation 
(load applied perpendicular to the rolling direction with the crack parallel to the rolling 
direction), stable crack growth will initiate in the interior. The crack will continue to 
grow stably at a higher rate in the interior than at the surface until the crack length has 
grown about 3mm at the surface. By this time, the crack front has developed a deeply 
tunneled shape that extends in a self-similar manner [17]. Shear lips did not form in this 
orientation. 

Thus, the initial stage of fracture in thin sheet aluminum appears to be dominated 
by variations in constraint, crack front shape, and CTOA. Furthermore, even in the 
steady state region, where the surface measurements of CTOA are constant, the crack 
front shape can be severely tunneled. These issues raise concerns over the validity of 
using surface measurements of CTOA as a fracture criterion in a two-dimensional 
analysis. 

The objective of this paper is to numerically examine the initial transient stage of 
fracture in thin- sheet materials as well as the effect of crack tunneling on the steady state 
region. A three-dimensional, elastic-plastic finite-element simulation of stable crack 
growth in thin sheet 2024-T3 aluminum alloy in the T-L orientation will be presented. 
The analyses considered both straight and tunneled crack fronts. Experimental data used 
in the numerical simulations were taken from Ref. 17. Comparisons are made between 
the critical crack-tip opening angle (CTOA), global constraint factor, and through- 
thickness stresses for the tunneled and straight crack fronts. 
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CRACK CONFIGURATION AND MATERIAL PROPERTIES 


A middle crack tension, M(T), specimen 76.2mm wide by 2.3mm thick (initial 
a/w=0.33 and a/B=5.56) was analyzed. The material was 2024-T3 aluminum alloy in the 
T-L orientation. The modulus of elasticity (E) was 71,400 MPa and Poisson's ratio (v) 

i 

was 0.3. In the analyses, a multi-linear representation of the uniaxial stress-strain curve, 
given in Table 1, was used. The stable tearing behavior was obtained from fracture tests 
where the applied stress (S) was recorded as a function of crack growth at the specimen 
surface (Aa s ) [17], as shown in Figure 1. An average of the applied stress against crack 
extension at the specimen surface data is given in Table 2. 

Deeply tunneled crack fronts were observed to develop during stable tearing tests 
of M(T) specimens [17]. The tunneled crack front shapes were obtained from fracture 
tests where the specimens were allowed to stably tear for a predetermined amount of 
crack growth at the specimen surface. Then, the specimens were fatigue cycled at a high 
stress ratio until failure. These changes in loading produced both macroscopic and 
microscopic changes in the fracture surface, creating a distinct boundary between the 
stable tearing and fatigue crack growth regions. Each boundary represents the crack front 
shape at a specific time during the stable tearing process. 

The shape of the crack front changed significantly during the early stages of stable 
tearing, as shown in Figure 2. Prior to stable tearing, the crack front shape was slightly 
tunneled due to the fatigue cycling used to generate a sharp crack-tip. During the initial 
stages of stable tearing, the crack grew faster in the interior. By the time the crack 
growth at the specimen surface reached 0.25mm, the growth in the interior was about 
1mm. The crack front shape stabilized after about 3mm of crack growth at the specimen 
surface. The stabilized shape was deeply tunneled,, also shown in Figure 2. Also, it was 
observed that the fracture surfaces were macroscopically perpendicular to the direction of 
loading (i.e., shear lips did not form) [17]. 

The crack front shapes given in Figure 2 were digitized and fit with polynomial 
expressions to describe the crack front shape as a function of through-thickness position. 
Additional crack front profiles were generated from the polynomial curve fits to describe 
the approximate crack front shape at each 0.25mm increment of crack growth at the 
specimen surface, as shown in Figure 3. 
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FINITE ELEMENT ANALYSES 


A three-dimensional, elastic-plastic finite element analysis code, ZIP3D [18] was 
used for the fracture simulation. The program uses 8-noded, hexahedral elements, the 
von Mises yield criterion, isotropic hardening, small-strain deformation theory, and 
Drucker's associated flow rule. The elastic-plastic analysis is based on the initial-stress 
method and the incremental theory of plasticity [18], The finite-element formulation and 
solution method are given in Ref. 19. 

Two different finite element models were created; one with straight, and one with 
tunneled crack front shapes. Symmetry conditions required that only one eighth of the 
M(T) specimens be modeled. Both models used five layers of elements through the half 
thickness (B/2=l. 143mm) and had 7398 nodes and 5650 elements. The thicknesses of 
the layers 1 through 5 (with the layer adjacent to the mid-plane denoted as layer 1) were 
0.125, 0.125, 0.125, 0.075, and 0.05 of B, respectively. The crack plane elements of the 
straight crack front mesh were 0.25mm in the y-direction and the x-direction, as shown in 
Figure 4. The crack plane elements of the tunneled crack front mesh were 0.25mm in the 
y-direction and in the x-direction at the surface, but the length of the x-direction elements 
in the interior were greater to reflect the extent of tunneling present. This is illustrated in 
Figure 3, where a view of the crack plane (x-z plane at y = 0) is shown. 

Normally, the nodes along the crack line and ahead of the crack- tip are fixed, 
while those behind are free. This analysis uses springs along the crack line to change 
boundary conditions associated with crack growth. The spring stiffness is set equal to 
zero for nodes behind the crack tip and assigned an extremely large value for nodes ahead 
of the crack tip. A nodal-release procedure [20] is employed for crack growth. Stable 
crack growth was simulated by allowing the crack to grow when the applied stress 
reached a predetermined level. These predetermined stress levels were obtained from the 
fracture test results reported in Ref. 17 and shown in Figure 1. In these tests, the crack 
length at the surface and the applied stress were measured during stable tearing. The 
curve fit shown in Figure 1 was evaluated to give the stress required for crack growth 
(critical stress) for each of the discrete crack lengths (measured at the surface) used in the 
finite element model. These values are listed in Table 2. In the finite element simulation, 
the entire crack front was allowed to advance one node when the applied stress reached 
the critical value corresponding to the current crack length at the surface. 
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RESULTS AND DISCUSSION 


Elastic-plastic finite element analyses were conducted to simulate fracture in t}ie 
M(T) specimens. Two models were examined; in the first, a straight through-thickness 
crack front was modeled and in the second, a tunneled through-thickness crack front was 
modeled. Loading was simulated with uniform tractions in the y-axes along the 
y=76.2mm surface of the models. The material stress-strain curve, applied stress-crack 
growth behavior at the specimen surface, and crack front profile shapes were obtained 
from experimental measurements [17]. The finite element analyses used this information 
to define the fracture behavior and generate the local crack-tip stresses, the global 
constraint factor, local crack-opening displacements, and through-thickness crack-tip 
opening angle (CTOA) values. 

Local Crack-Tip Stresses 

The stresses ahead of the crack tip in the direction of loading (o yy ) were obtained 
from the finite element analyses of the straight and tunneled crack fronts. The reported 
cj yy stresses are at the applied stress level required for crack growth, but before the crack 
front was allowed to grow. Due to the singularity at the crack tip and the approximations 
inherent in the small strain finite element analysis, the stresses in the crack-tip element 
were ignored. However, the stresses in the other elements should be reasonably accurate. 

The Oyy stresses, normalized by the flow stress (Go), on the crack plane ahead of 
the crack-tip in the straight crack front analysis are shown in Figures 5a and 5b for the 
surface and mid-plane elements, respectively. Prior to crack growth, the trends of the Oyy 
stresses agree with the large deformation theory finite element results of Horn and 
McMeeking [14]. Their analysis found that, for a blunt notch in a thin sheet, the Oyy 
stresses at the mid-plane were greater than those at the surface with normalized peak 
stresses (Oyy/Oo) of about 1.2 at the surface and 2.0 at the mid-plane. The present results 
exhibited similar behavior with normalized peak stresses of 1.3 at the surface and 2.0 at 
the mid-plane. 

The Oyy stresses, normalized by the flow stress (Oo), on the crack plane ahead of 
the crack tip in the tunneled crack front analysis are shown in Figures 6a and 6b for the 
surface and mid-plane elements, respectively. For the tunneled crack front, the trends of 
the normalized Oyy stresses at the mid-plane as a function of distance from the crack tip 
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agree with the straight crack results for the mid-plane. However, the surface results are 
different for the two crack front shapes. After crack growth, the normalized stresses 
ahead of the crack-tip in the tunneled crack front model decrease from a value around 2, 

reach a minimum of about 1 and increase to a second peak before decreasing once again. 

< 

The differences in the normalized a yy stresses along the crack plane in the 
tunneled and straight crack front finite element analyses are further illustrated in the 
contour plots of Figure 7. These plots show lines of constant o yy /oo stress in the steady- 
state region of crack growth (just prior to the 10th increment of crack growth) for an area 
extending from the crack front to about 3mm ahead of the crack front. Also 
superimposed on these plots are the finite element sizes along the crack plane. The high 
Cyy/do stress gradients occur over a minimum of 5 elements, which, though not sufficient 
to capture the magnitude stresses at the crack front, should be sufficient to capture the 
behavioral trends and the stress magnitude several elements ahead of the crack front. The 
finite element results indicate that the wedge formed at the surface of the tunneled crack 
front significantly increases the near surface stresses. The stresses along the mid-plane of 
the tunneled crack front are also greater than those at a similar position in the straight 
crack front analysis. 

Global Constraint Factor 

The global constraint factor (a g ) was defined by Newman, Bigelow, and 
Shivakumar [15] to be the average normal-stress-to-flow-stress ratio for elements in the 
plastic zone. 



where, A m is the projected area on the uncracked ligament of a yielded element m. At is 
the total projected area for all elements that have yielded, and M is the total number of 
yielded elements. Newman et al [16] found the global constraint factor to be nearly 
independent of specimen type for crack length to thickness (a/B) and uncracked ligament 
to thickness ratios (b/B) greater than 4. An equation expressing the global constraint as a 
function of the stress-intensity factor was obtained from elastic-plastic finite element 
analyses of stationary, straight through-thickness cracks in an elastic-perfectly plastic 
material [16] as 
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t 

Using Eqn. 1, the global constraint factor for the straight and tunneled crack 
analyses are plotted in Figure 8 as a function of the normalized stress-intensity factor. 
For the tunneled crack, the stress-intensity factor was calculated, using Eqn. 4, neglecting 
the crack front curvature and using the surface crack length measurement (Note: that by 
using the surface crack length, the average stress intensity factor of the tunneled crack 
would be about 10% less than that obtained using the crack length at maximum depth. 
This surface crack representation was used throughout the paper to be consistent with 
experimental measurements, where only the surface value is known.). Initially, a g for the 
tunneled crack is slightly greater than that of the straight crack for the same normalized 
stress-intensity factor. As the stress-intensity factor increases, a g levels off at 1.3 for 
both the crack front shapes. This value is higher than that calculated from Eqn. 2, 
however, the present analysis simulates a material that strain hardens rather than an 
elastic-perfectly plastic material. As the crack begins to grow, the global constraint factor 
drops until the maximum stress is reached. A slight increase in a g was observed after the 
peak stress was reached. 

Crack Surface Displacements and CTO A 

During the finite element simulation, the displacements of the nodes behind the 
crack tip on the crack plane were recorded at several stress levels and prior to each 
increment of crack growth. The displacements in the direction of loading (v) were used 
to calculate the crack-tip opening angle (CTOA). CTOA was calculated in the same 
manner as Dawicke and Sutton [9] using surface measurements. CTOA is defined to be 
the angle (\jr c ) made by a point on the upper crack surface, the crack tip, and a point on 
the lower crack surface directly below the one on the upper surface. The locations of the 
points on the upper and lower crack surfaces were 1mm behind the crack- tip if the 
amount of stable tearing exceeded 1mm, or at the node ahead of the original crack front if 
the total amount of stable tearing was less than 1mm. This is illustrated in Figure 9 for 
displacements obtained from the finite element analysis. This definition was used to 
prevent anomalous CTOA values that would result from the blunted region formed at the 
original crack-tip. In terms of the finite element results, \|/ c is determined from the v- 
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displacement at a distance 1mm behind the crack- tip (or at the node immediately ahead of 
the original crack tip if Aa < 1mm) and the x-distance from the crack- tip to the point of v- 
displacement measurement. 

Vc = 2tair'(|) 1 (3) 

The CTOA was calculated just prior to each increment of crack growth (except 
the first) for both the tunneled crack front and the straight crack front models, as shown in 
Figures 10a and 10b, respectively. CTOA values were calculated at the surface 
(z/B=0.5), the mid-plane (z/B=0), and at different through-thickness locations 
(z/B=0.125, 0.25, 0.375, and 0.45). The tunneled crack front simulation exhibited a 
considerable variation in CTOA through-the-thickness during the initial 3mm of crack 
growth at the specimen's surface. The surface CTOA values agreed with the 
experimental surface measurements, with an initial angle of 10° at a crack growth at the 
specimen's surface of 0.5mm. The calculated CTOA values decreased through-the- 
thickness to a value of about 4° at the mid-plane. After 3mm of crack growth at the 
specimen's surface, all through-thickness calculations of CTOA approached the average 
CTOA (4.7°) that was experimentally measured at the surface [17]. The straight crack 
simulation showed very little variation in CTOA through-the-thickness. Initially, the 
calculated CTOA value was between 5° and 6° (with the lower value at the surface). As 
the simulation grew the crack, the calculated CTOA values decreased slightly to between 
4° and 5°. 

To further illustrate the through-thickness variation in CTOA for the tunneled 
crack front, the crack opening displacement was plotted as a function of crack length and 
position through-the-thickness, as shown in Figures 11a and lib for crack growth at the 
specimen surface of 0.75mm and 4.75mm, respectively. The crack opening 
displacements shown in Figure 1 1 are just prior to the next increment of crack growth. In 
Figure 11a, the crack has grown three nodes past the original crack- tip. The large initial 
tunneling in the interior is evident, as is the severe blunting of the original crack front and 
the variation in CTOA through-the-thickness. The first two elements along the mid-plane 
are significantly larger (in the crack growth direction) than those at the surface due to the 
severe tunneling during the early stages of stable tearing. However, after the second 
increment of crack growth, the crack growth at the mid-plane becomes closer to that at 
the surface, yet the differences in CTOA remain large (see Figure 10a). In Figure 11b, 
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the crack has reached a "steady-state" condition, where the crack front profile (imposed 
condition) and the CTOA (calculated result) are constant through-the-thickness. 

A convergence study was not performed to assess the effect of element size on the 
through-the-thickness CTOA behavior. However, tfye element size used in this study was 
small enough to provide accurate predictions of stable tearing and local crack opening 
displacements in earlier two-dimensional predictions [8-9, 11-13]. The CTOA was 
determined four elements away from the crack front (except for the second and third 
increments of crack growth where two and three elements were used, respectively). 
Additional layers through-the-thickness and smaller elements would give more accurate 
results, but the trends shown in this analysis would probably not change. 

The tunneled crack front finite element simulation indicates that CTOA may be a 
fracture criterion for thin-sheet materials that can be incorporated into a two-dimensional 
analysis. The three-dimensional analyses indicate that after steady-state conditions have 
been achieved, the CTOA measured at the surface is representative of what is happening 
in the interior. In the transient region the surface CTOA measurements are higher than in 
the interior and the average CTOA through-the-thickness is close to the steady-state 
CTOA. In a two-dimensional fracture simulation, using the constant steady-state CTOA 
during the transient region may be an adequate approximation, especially if the amount of 
stable tearing is large compared to the transient region. Earlier two-dimensional, elastic- 
plastic finite element analyses, using the measured steady-state CTOA value as a fracture 
criterion, have successfully predicted fracture behavior of single and multiple cracks in 
2.3mm thick by 300mm wide sheets of 2024-T3 [21]. 


CONCLUSIONS 

A small strain theory, three-dimensional elastic-plastic finite element analysis was 
used to simulate fracture in thin sheet 2024-T3 aluminum alloy in the T-L orientation. 
Both straight and tunneled crack fronts were modeled. The tunneled crack front shapes 
were obtained from experimental measurements. The stable crack growth behavior was 
measured at the specimen surface as a function of applied stress. The fracture simulation 
modeled the crack tunneling and extension as a function of applied stress. The results 
indicate that: 
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(1) The global constraint factor, a g , dropped at the onset of stable crack growth. 
After peak stress was achieved, a g began to increase. The effect of crack front 
shape on a g was small. 

(2) The crack front shape did influence botjh the local constraint and through-the- 
thickness CTOA behavior. Failure to account for these crack shapes in any three- 
dimensional simulations or predictions may give inaccurate results. 

(a) Local constraint differences: Prior to crack growth, the local a yy stresses 
ahead of the crack-tip were higher at the mid-plane and lower at the surface 
in the tunneled crack model than they were in the straight crack model. In 
the straight crack model, both the surface and mid-plane <T yy stresses 
increased slightly with crack growth. In the tunneled crack model, the local 
stresses dropped during the initial crack growth. During the rest of the initial 
transient region (Aa s < 3mm), the mid-plane stresses increased sharply and 
the surface stress only slightly increased. After the initial transient region, 
the mid-plane stresses in the tunneled crack model were greater than and the 
surface stresses were less than those of the straight crack model. 

(b) CTO A differences: Results from the tunneled crack front model indicated 
that during the initial transient region (Aa s < 3mm) the CTOA at the mid- 
plane was less than 4° and increased with crack growth. The initial CTOA at 
the surface was about 10° and decreased with crack growth, which agreed 
with experimental surface CTOA measurements. At other through-thickness 
locations, the CTOA values fell between the surface and mid-plane values. 
After the transient region, the experimentally measured through-the-thickness 
CTOA was nearly constant, at a value of 4.7°. For the straight crack front 
model, the CTOA was constant thought-the-thiekness and the initial surface 
CTOA values did not agree with the experimental measurements. 

(3) The three-dimensional analysis appears to confirm the potential of CTOA as a 
two-dimensional fracture criterion. In the initial transient region, the high angles 
measured at the surface do not reflect the behavior in the interior and an average 
value might be close to the steady state surface measurements. In the steady state 
crack growth region, the CTOA appears to be constant though-the-thickness. 
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Table 1 

Multi-linear representation of the uniaxial stress-strain curve for 2024-T3 aluminum alloy 

in the T-L orientation [17] 
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Table 2 

Stable crack growth data for 2024-T3 aluminum alloy in the T-L orientation [17] 
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Experimental measurements of stress against crack growth at the surface 
and the multi-linear curve fit for 2024-T3 aluminum alloy in the T-L 
orientation 
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Figure 2 Crack front profiles during stable tearing for 2024-T3 in the T-L 
orientation [17]. 
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Figure 3 Crack front profiles at even increments of crack growth at the specimen 
surface, generated from polynomial curve fits to the measurements given 
in Figure 2. 
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Finite element mesh model of l/8th of the M(T) specimen 
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a. Surface (z/B=0.5) 
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b. Mid-plane (z/B=0) 

Figure 5 Normalized a yy stresses at the second element ahead of the crack-tip from 
the elastic -plastic finite element analysis of the M(T) specimen with a 
straight crack front. 


19 




Distance Ahead of Crack Tip 

(mm) 


a. Surface (z/B=0.5) 



Distance Ahead of Crack Tip 

(mm) 
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Figure 6 Normalized o yy stresses (after the first element ahead of the crack-tip) 
from the elastic-plastic finite element analysis of the M(T) specimen with 
a tunneled crack front. 


20 



a 3.5 
b 3.25 
c 3.0 
d 2.75 
e 2.5 
f 2.25 
g 2-0 
h 1.75 
i 1.5 
j 125 


z 

L. 


a\\c' 


a ^ 
yy cor 
/ C7_ li 


lmm" 


contourj 

lines 


a. Curved crack finite element analysis 
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b. Straight crack finite element analysis 

Contour plots of Oyy/o 0 ahead of the crack front before the 10th 
increments of crack growth and the mesh pattern along the crack plane. 
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Figure 8 Global constraint factor, ag, as a function of the normalized applied stress- 
intensity factor. 
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a. 0.75mm of crack growth, \|/ c measured 0.5mm behind crack- tip 
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b. 3mm of crack growth, \|/ c measured 1mm behind crack-tip 

Figure 9 Crack opening displacements at the surface for a tunneled crack front 
indicating the location of CTOA measurement. 


23 




(mm) 


a. finite element results for tunneled crack front 



Aa s 


(mm) 

b. finite element results for straight crack front’ 

Figure 10 Experimental surface crack-tip opening angle measurements [17] and 
finite element results for a 2024-T3 aluminum alloy in the T-L orientation. 
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a. Aa s = 0.75mm 


Figure 11 



(mm) 


b. Aa s = 4.75mm 


Through-the-thickness crack opening displacement profiles for a tunneled 
crack front at two stages of crack growth. 
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